Astronomy &: Astrophysics manuscript no. h3741 


February 2, 2008 


(DOI: will be inserted by hand later) 





Spectral modeling of the three spectral states of the Galactic 

microquasar GRS 1915+105 

Malgorzata A. Sobolewska and Piotr T. Zycki 

Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 
e-mail: (malsob, ptz)@caink. edu.pl 

Received accepted ... 

Abstract. 

We have analyzed Rossi X-ray Timing Explorer spectral data of the Calactic micro-quasar CRS 1915-f 105 in 
its various spectral states, as defined by Belloni et al. (2000). In states A and B the spectra are dominated by 
a strong soft thermal component, accompanied by a weak harder tail. The soft component is rather complex 
and cannot be described as a simple accretion disk emission. Relativistic effects in Kerr metric contribute to the 
complexity of the soft component but are not sufficient to fully account for it. As found previously, state C spectra 
are dominated by a Comptonized component, with small contribution from disk photons. The X-ray reprocessed 
component is highly significant in those spectra and, in contrast to the usual hard state spectra from accreting 
black holes, it is highly ionized. 
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1. Introduction 

Ever since its discovery by Castro-Tirado et al. H1994|l . 
. the transient Galactic source GRS 1915-1-105 has been a 
' subject of intense observational and theoretical work. The 
, source shows a number of remarkable features, most spec- 
' tacular of which are ejections of relativistic jets with ap- 
, parent super-luminal motion (Mirabel & Rodriguez ll994|l . 
■ It displays a large variety of behaviour in the X-ray band 
, (dramatic flux and spectral variability on time scales from 
' miliseconds to months) , most likely reflecting violent phe- 
nomena occurring in the inner accretion flow. 

Large amount of data collected by various X-ray ob- 
servatories indicate that the source stays permanently in 
a soft spectral state (see, e.g., Nowak|1995 for a charac- 
teristics of the various spectral/ timing states). That is, its 
spectral energy distribution is dominated by the 1-10 keV 
band, even though the ratio of soft thermal (temperature 
~ 1 keV) emission to harder Comptonized continuum does 
vary substantially. The source never seems to have been 
in the hard spectral state, which is more typical for black 
hole binaries with energy distribution peaking at ^ 100 
keV (see Pout anen 1 1 9 9 81 for review). 

A characteristic feature of GRS 1915-1-105 is its strong 
variability, of both flux (count rate) and energy spectra. 
Among the variety of light curves and color-color dia- 
grams certain patterns can be found. These were stud- 
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icd in detail by Belloni et al. (2000; hereafter IB00|l . using 
RXTE data, who attempted to classify the behaviour of 
the source. They divided the lightcurves into 12 classes 
based on the color-color diagrams and count rates, and 
found that the lightcurves from each class can be decom- 
posed into only 3 basic spectral states (A, B, C in the 
nomenclature of IBOO|l . The strongly variable lightcurves 
of the source could be explained in terms of rapid switches 
between these three states. 

There have been previous attempts to model the X-ray 
spectra of GRS 1915+105 in the B and C spectral states. 
Vilhu et al. l)^niIT|l analyzed RXTE data (PCA/HEXTE 
instruments, in the range of 2-60 keV) and concluded that 
spectra of GRS 1915-1-105 can be described with the model 
consisting of a disc blackbody and a Comptonized compo- 
nent with Compton reflection. Zdziarski et al. ( 2 001|l pre- 
sented an analysis of simultaneous observations of GRS 
1915+105 in states B and C by RXTE and CGi?0/OSSE 
instruments. They found that the broad-band spectra re- 
quired non-thermal Comptonization of soft photons, since 
the data do not show any break up to at least 600 keV. Low 
temperature, thermal Comptonization was also needed to 
correctly model (together with blackbody component) the 
soft X-ray part [E < 10 keV), thus the best final model 
invoked a hybrid plasma (Coppi ,1999p . Vadavale et al. 
1)200 l|l focused on study of the source during the steady 
low-hard state (class x, state C), which they further subdi- 
vided into radio- loud (xrl) and radio-quiet (xrq) ones. 
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From their analysis followed that the spectrum of GRS 
1915+105 in state C can be described as a superposition of 
disk blackbody and Comptonized component, and an ad- 
ditional power low component in the case of the xrl state. 
This additional power low was interpreted as an emission 
from the base of a jet. Analysis of spectral features due 
to iron emission/ absorption in the 5-9 keV band was per- 
formed by Martocchia et al. (|2)02) using BeppoSAX data 
of the source in state C, according to their interpretation. 
They found that the Fe Kq, line is broad and it has the 
characteristic asymmetric profile indicating its origin from 
a rotating accretion disk. 

In all above mentioned analyses, disk blackbody or 
blackbody component is used to (fully or partially) de- 
scribe the soft emission from the disk. The diskbh model 
reproduces the shape of the disk spectrum quite accurately 
unless the relativistic corrections are extreme. However, 
the numerical values of the parameters do not easily trans- 
late into physical parameters (Merloni et al. l2(JU(Jfl . These 
can be estimated from a relativistic model, which includes 
radiative transfer effects, at least approximating them by 
introducing the color temperature correction (Shimura & 
TakaharaEMSI)- 

Therefore, in this paper we analyse data from all three 
spectral states of GRS 1915+105 , applying the same set 
of models to all the data. The models include one for emis- 
sion from an accretion disk around a maximally rotating 
Kerr black hole. 

2. Data selection and reduction 

The log of PCA/HEXTE observations used for our anal- 
ysis is given in Table ^ The data were extracted from 
pubhc archive at HEASARC/GSFC and reduced with the 
FTOOLS software package (version 5.0). 

Based on BOO we attempted to extract spectra corre- 
sponding to 'pure' spectral states A, B and C. We have two 
observations from class x (pure state C), one observation 
from class A from which we extract a period corresponding 
to state B based on the X-ray color-color diagrams and 
lightcurves appearance fsee lBOOl for details), one observa- 
tion from class 9 which provides state A spectra, and two 
observations from class (3 - further examples of state A, 
B, and C spectra. 

For the spectral analysis, PGA Standard 2 mode data 
(PCU 0-2, top layer only) and HEXTE Archive Mode 
(both clusters) configuration are used. The energy range 
is 3-20 keV and 20-100 keV, correspondingly. We assume 
systematic errors in the PGA data at the level of 1% due 
to the uncertainty of response matrices. Standard back- 
ground subtraction and dead-time correction procedures 
are applied. 

The value of the hydrogen column density, Nh, is al- 
lowed to vary during the fit. We assume the source in- 
clination is J = 70°, distance d = 12.5 kpc (Mirabel & 
Rodriguez '1994) and central black hole mass M — 14Mq 
(Greiner et al. i2001|l . We allow for free relative normaliza- 



Table 1. The log of observations. The letters I, K 
stand for observations ID 10408-01, 20402-01, respectively 
(|B001I1 ■ The third column indicates which state spectra 
were extracted from the data. 



Observation ID 


Class 


State 


K-44-00 


/3 


A, B, C 


K-45-03 


/3 


A, B, C 


K-45-02 


6 


A, B 


1-38-00 


A 


B 


K-05-00 


x' 


C 


1-42-00 




C 



tion of HEXTE data with respect to PGA data. Analysis 
was performed in XSPEC ver. 11 fArnaud ll996|l . 

3. IVIodels 

The disk soft thermal emission will be described by the 
simple disk blackbody model (Mitsuda et al. 119841 im- 
plemented as diskbh) in XSPEG, or a relativistic accre- 
tion disk model for Kerr geometry (Ebisawa et al. 12001(1 . 
The KerrD model is computed for the maximally rotating 
black hole (a = 0.998; Thorne [Tn71|l . It assumes blackbody 
spectra at each radius with the color temperature correc- 
tion constant with radius and convolves them with the 
transfer function describing photon propagation, as com- 
puted by Laor H1991|l . We fix the black hole mass, source 
distance and accretion disk inclination at the values indi- 
cated above, and allow only the inner radius of the disk 
(in units of i?g = GM /c^) and the mass accretion rate, M 
to be free parameters. In particular, the model does not 
have an additional free normalization. 

The Gomptonized component will be described by the 
thComp model, based on a solution of the Kompaneets 
equation (Zdziarski et al. ll996|l . or by the compPS model 
of Poutanen & Svensson fl996). The latter does not as- 
sume the diffusion approximation but solves for each scat- 
tering order independently. The compPS model can be 
used with purely thermal, purely non-thermal or a hybrid 
electron energy distribution (Goppi 1999). The seed pho- 
tons will be assumed to have a disk blackbody spectrum 
and the comptonizing cloud is assumed to be a sphere. The 
code also consistently computes the un-scattered fraction 
of seed photons. 

For the X-ray reprocessed component we will use the 
relrepr model (Zycki et al. 12001(1 . combining the Gompton- 
reflected continuum (Magdziarz & Zdziarski 119951 Done 
et al. T992") with the Fe K^ line computations of Zycki & 
Gzerny (1994). The model computes self-consistently the 
strength of the line for a given primary continuum spec- 
trum, amplitude of the Gompton-reflected component, R, 
ionization parameter, ^ — AttFx /n, and Fe abundance, 
and it also allows for relativistic effects to be applied to 
both components (using prescription from Fabian et al. 
ri989)l . The amplitude R is defined as n/{2n), where f2 
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Table 2. Values of , {xD for various models fitted to the data. 



State Obs. Id (1) (2) (3) (4) 





K-44-00 


liM(i.i4) - 






c 


K-45-03 


W (0-83) - 








K-05-00 


W (0.87) - 








1-42-00 


lt2(i.04) - 







K-44-00 4.8 (0.70) %S (o.68) ^ (0.66) 

A K-45-03 3.4 (1.23) ^ (1.19) ^^ (1.29) 

K-45-02 3.4 ^ (0.83) ^ (0.83) ^ (0.96) 



K-44-00 2.7 %a (0.79) %i (0.70) ^ (0.65) 

B K-45-03 1.6 %i (1.00) (0.97) ^ (0.93) 

1-38-00 2.3 ^ (0.93) ^ (0.89) (0.84) 

Models (2), (3) and (4) fitted to state A spectra contain an additional gaussian line 
Fitted models: 

(1) compPS - comptonized disk blackbody 

(2) thComp +compPS - two comptonized disk blackbody spectra 

(3) hbody +compPS - comptonized disk blackbody with an additional blackbody 

(4) KerrD+compPS - Kerr metric disk spectrum and a comptonized disk blackbody 
In all fits the compPS model includes the reprocessed component 



is the solid angle subtended by the reprocessor from the 
X-ray source. 

4. Results 

4.1. State B spectra 

We begin our analysis with the brightest state, B. First, 
we try modeling soft part of the spectrum with disk black- 
body, and the hard one with non-thermal Comptonization 
(we apply compPS model to the data) . The latter is moti- 
vated by the reports, where the source spectrum in state 
B and C was found to extend to > 600 keV without a 
break or a cutoff (Zdziarski et al. I2(JU1|I . We assume that 
the seed photons for Comptonization have the same disk 
blackbody shape and temperature. The reprocessed com- 
ponent is also included in the model. This model gives 
unacceptable fits, xt ~ 2, with residuals up to 10 per cent 
(Model 1 in Table|2I). Good description of data is provided 
by models including a more complex shape of the soft 
component. Motivated by earlier results (see Zycki et al. 
121)011 Wilson & Done EUrm and references therein) we test 
two phenomcnological descriptions of the soft component: 

(i) with an additional Comptonized disk blackbody and 

(ii) an additional blackbody. (Note that a certain fraction 
of the seed disk blackbody photons always contributes to 
the soft component due to modeling the hard part of the 
spectrum with compPS .) Both these complex models give 
significant improvement compared to Model 1, with best 
fits of comparable quality, as shown in Table El (Models 2 
and 3). 

A possible physical reason for the complexity may be 
the relativistic effects expected from an inclined accretion 



disk. Therefore we also try (iii) the KerrD model for the 
soft component (with non-thermal compPS for the hard 
tail, as before). The quality of the fits is similar for all 
the three models (Model 4 in Table E)). The relativistic 
effects are quite pronounced at the high inclination of the 
source, i « 70° (see Ebisawa et al. 120011 Gierlinski et al. 
12001(1 . and the Doppler-boosted disk emission fits well the 
observed high energy cutoff of the soft component. In the 
best fit for K-44-00 te inner disk radius i?in = 3.80^2:56 -Rg> 
for i — 70° and d — 12.5 kpc. With the more recently 
determined i = 66° and c? = 11 kpc (Fender et al. I1999|l 
the inner radius is i?in = 2.25t.i°i- In both cases Rin is 
consistent with the last stable orbit at 1.24 Rg, and a large 
inner hole in the disk is ruled out. A good fit jv — 
46.3/72) can be obtained with the disk emission model in 
Schwarzschild metric, however a small mass of the black 
hole, M w 4M0 is required in this case, in order to account 
for the high observed disk temperature. Thus, if the soft 
component is modeled as a relativistic disk emission, a 
Kerr (rotating) black hole is implied by our result. 

We note once again that some seed disk blackbody 
photons do contribute to the soft component (Fig. ^) 
since the compPS model computes self-consistently the 
un-scattered fraction of the seed photon flux. Without 
this component the fit is worse by Ax^ = 21 (for K-44- 
00) , therefore we stress that the relativistic effects may be 
not the only factor contributing to the complexity of the 
soft component. However, given the high inclination of the 
accretion disc, their a priori importance is obvious. The 
mass accretion rate is 7.5 x 10^* gs^^ for K-44-00 data 
(Table El and up to 2.6 x lO^^gs^^ for other data). We 
note that both M and -Rin would be different for different 
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value of the black hole angular momentum (see discussion 
in Gierlinski et al. 120011 recall that a = 0.998 is assumed 
in our paper), thus our result can only be regarded as a 
qualitative demonstration that relativistic effects can par- 
tially account for apparently complex soft component in 
this spectral state. 

The seed photon temperature for the non-thermal 
Comptonization was a free parameter in these fits, but 
the best fit values are similar to the obtained maximum 
temperatures of the Kerr disk emission (« 1.2 keV), 
as expected. The reprocessed component from the non- 
thermal primary is present in the data, with typical am- 
plitude R — 0.7-1.0 (TableEJ. Its presence is however not 
highly significant, since removing it we obtain fits worse by 
Ax^ — 5.4-8.4. The best fit value of ionization parameter 
is high, ^ ~ 10*, as expected from a hot accretion disk, al- 
though the formal errors on ^ are large. No further broad- 
ening of the reprocessed component is formally required 
by the data, although the data are consistent with rela- 
tivistic smearing from a disk extending down to ~ 20i?g. 
This is larger than the inner radius of the emitting disk, 
indicating that the innermost disk region is too hot to 
produce any significant spectral features. A part of the 
uncertainty in determining i?in comes from the fact that 
Compton scattering in highly ionized plasma contributes 
to smearing of the spectral features, making quantitative 
determinations of Rm rather uncertain (Ross et al. ll999|) . 

4.2. State A spectra 

State A is very similar to state B, although somewhat 
less luminous (jBOOjl . Presence of the non-thermal tail was 
not explicitly demonstrated in this state, but given its 
similarity to state B, we assume the hard tail is indeed 
given by a non-thermal Comptonization and include it 
in all models. Initial fits with the non-thermal compPS 
model give bad fits {xl ~ 3-5). Again, a complex shape of 
soft component helps to obtain good fits, with comparable 
quality of fits for the various models (Table 0). The overall 
spectra are dominated by the strong soft component, with 
the harder tail appearing only above ~ 20 keV (Fig. ^ 
Table El. 

The model which includes the Kerr metric disk com- 
ponent gives comparable fit quality to the other models 
with complex soft component. Moreover, the unscattcred 
photons contribute significantly to the soft component, 
broadening it compared to what would be obtained with 
pure relativistic disk emission (Fig. Removing the un- 
scatered photons from the model gives clearly unaccept- 
able fits (xi ^ 3), i.e. the Kerr disk component with pure 
non-thermal Comptonized component cannot describe the 
data. In one case (K-45-03) the best overall fit is obtained 
with the blackbody model while in another case (K-45- 
02) the Comptonized disk blackbody and the blackbody 
models provide the same quality of fits to the data. The 
temperature of the additional blackbody is rather high, 
~ 2 keV. It appears that this component is used to adjust 



the curvature of the Comptonized spectrum rather than 
represents a real spectral component. Clearly, relativistic 
effects from a Kerr metric may be considered as contribut- 
ing to the complexity of the soft component in the case of 
K-44-00 observation, but are not sufficient to describe it 
completely. 

Other aspects of the best fit models are similar to state 
B spectra: the reprocessed component is present in the 
data, it is highly ionized and further smeared (again, with 
rather large errors). An additional narrow gaussian line 
at 6.4 keV was also included in the model, since the fit 
residuals clearly indicated the presence of such a line. Its 
equivalent width is 30 eV. 

Fits of the KerrD model with the non-thermal tail give 
M systematically lower than those for state B, while i?in 
is the same within the fit unccr tainty (Table OJ. 

4.3. State C spectra 

State C is generally rather harder than either A and B, 
and the disk component is less pronounced ljBOO|l . The 
non-thermal compPS model gives adequate fits to the data 
(with our assumed systematic errors of 1%). The xt val- 
ues for K-44-00, K-45-03, K-05-00, T42-00 observations 
are xl = 1-14, 0.83, 0.87 and 1.04, respectively (Table|21 
model 1; Fig. HI). 

The disk component can be separated from the total 
Comptonization component as computed by compPS by 
setting the un-scattered fraction of seed photons to 0, and 
modeling the soft photons by a separate model. We per- 
form such a procedure for the K-44-00 and K-45-03 data 
sets, describing the soft component by the KerrD model. 
Although it does not lead to satisfactory fits, it allows to 
approximate values of physical disk parameters, M and 
i?in. These are A/ « 2 x lO^^gs^^ and R-.n « 15 Rg. The 
mass accretion rate is rather high, close to the upper val- 
ues for state B. The large inner radius means that the 
maximum color temperature is « 0.8 keV and most of the 
corresponding luminosity is not seen in the X-ray band. 

In all data sets the reprocessed component is always 
significantly present. Its amplitude is smaller than 1, it is 
highly ionized (^ ~ 10"*, but with large errors). Additional 
broadening of the spectral features is also required. If this 
is interpreted as relativistic, the inner radius of the reflect- 
ing disk is - 25-100 i?g (TableEJ. 

5. Discussion 

Our analysis of RXTE data of the micro-quasar GRS 
1915-1-105 reveals all the spectral features typically found 
in high/soft states of accreting black holes (see Gierlinski 
et al. nM?l Zycki et al. [2001: Zdziarski et al. IM7T]l . 
The hard component is well fit by a non-thermal 
Comptonization, the soft component is usually complex, 
and the reprocessed component is present, highly ionized 
and further broadened and smeared. 

Clear complex soft component is seen in the high lu- 
minosity states A and B. A number of descriptions of the 
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Fig. 1. KerrD model fits to spectra from observation K-44-00. Long-dashed curves are the comptonized component 
compPS with the transmitted seed photons, Short-dashed curves are the Kerr disk spectrum model, and the dash-dot 
curves are the reprocessed component. Complexity of the soft component in states A and B is apparent and this cannot 
be accounted for by the Kerr disk model alone. There is only a weak disk thermal component in state C spectrum, 
but modeling it with the KerrD model suggests that the disk M may actually be comparable or even higher than 
in states A and B. However, the inner disk radius is larger than the last stable orbit, so most of the accretion power 
would not appear as thermal X~ray emission. 



Table 3. Fit parameters and values for observation K-44-00 in states A, B, and C. Model consists of KerrD and 
non-thermal compPS (A and B), or single compPS (C) 



Parameter 


State B 


State A 


State C 


riH [xlO^^ cm~^] 


"•■^-0.2 




r, (i+0.7 


M [xlO^* gs"^] 








-Rin [-Rg] 


3.80t^;f, 


J^-^' -0.23 




fcTin (keV) 




l-16tH^3 




Fe 






O-'-0.5 


T 


U.UO_Q 05 


nfi+°°* 

U-UD_o.04 


0.441°:- 


R— ^ 

zn 


1 (fixed) 


0.70tg°53 


0.30i°;?° 




2200t^^r 


1.10tJ°o2XlO* 


(5200+23000 

OZUU_gooo 


-Rin.rcfl [-Rg] 


1651^45 


i37tr3i 






0.65 


0.66 


1.14 C-W) 



Parameters: 

M - mass accretion rate in the KerrD model 

-Rin - inner radius of the disk in the KerrD model in -Rg — GM/(? units 

fcTin - seed photons temperature for Comptonization in keV 

Fe - index of a power law electron energy distribution 

r - Thomson depth of the Comptonizing cloud 

Q./2iT - solid angle of the reprocessor 

^ = 4:TTFx/ne - ionization parameter of the reprocessor 

Rrafi - inner radius of the reprocessing disk 



component are possible: in addition to the usual disk ther- 
mal emission, there may be contribution from hotter spots 
on the disk, heated by strong illumination by compact 
magnetic flares. This would produce an additional black- 
body emission, such as in our model 3. Alternatively, the 
hot upper layers of the heated accretion disk may addi- 
tionally comptonize the soft photons from below, produc- 
ing the second, low temperature comptonized component 
(see more detailed discussion in, e.g., Zycki et al. 2001). 
However, the relativistic distortions expected from such a 
highly inclined accretion disk {i « 70°) are significant and 



they cannot be ignored. In particular Doppler blueshift 
and boosting make the observed temperature higher than 
the true color disk temperature. However, even assuming 
a maximally rotating black hole, we need to include a con- 
tribution from additional soft photons, in order to explain 
the spectrum in both state A and B. Therefore, quantita- 
tive determinations of the black hole angular momentum 
are not possible without constructing a more detailed ge- 
ometrical model explaining the origin of all components. 
However, relativistic effects alone in Schwarzschild metric 
(a = 0) would not be able to account for the high apparent 
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peak energy of the soft component, if the mass of 14 Mq is 
adopted. In the interpretation of Zdziarski et al. ( 2001) the 
high apparent peak energy is due to Compton upscatter- 
ing by a relatively cool, optically thick plasma (see fig. 3b 
in their paper), but clearly the Comptonization parame- 
ters will be changed if relativistic effects are taken into 
account in this description. 

In state C we generally see strongly comptonized emis- 
sion with small contribution of soft seed photons. The pa- 
rameters of the soft component correspond to larger inner 
disk radius than in states A and B, confirming qualita- 
tively previous studies of Belloni et al. H1H97(I . We note 
however that quantitative estimates of e.g. the inner disk 
radius require proper decomposition of spectra into the 
disk and Comptonized emission. Models used in many pre- 
vious papers, i.e. diskbb and a power law are clearly in- 
adequate for such a task. In particular a power law is not 
a proper description of a Comptonization process, if the 
seed photons are visible in the considered energy band, 
as is clearly the case in GRS 1915-1-105. Our modelling 
suggests inner radii « 15 -Rg. Since in a Kerr metric the 
dissipation of gravitational energy is strongly concentrated 
towards the center, therefore the total accretion power is 
implied to be much larger than that directly observed as 
an X-ray emission. Quantitative estimates would require 
knowledge of the black hole spin, but even for a mild value, 
a — 0.5, the fraction of energy dissipated within 15i?g is 
« 80%. Total accretion power implied would then not be 
different in state C compared to states A and B, but most 
of that power would be used to e.g. cause an ejection of 
plasma. 
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